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Linear programming models are useful tools in energy systems analysis 


Throug' 
oug terfuel substitutability, changes in Price structure affect pat- 


ር 
erns of resource utilization. However, total energy demand and supply c 
on 


Strai 
nts are deterministic. In this paper, several methods of representing 


supply and demand elasticities are considered. It is shown that when th 
e 


demand (o ! i 
(or supply) curves are Plecewise linear, the equilibrium demand ۲ 


wh ۱ 
ere the demand curve prices equal the marginal costs, can be found by sol 
b ۷ 


a d 
quadratic programming problem. In some cases these equilibrium points can 


b 
e determined through Simple modification of the linear programming algorith 


ing 


— , <۰ تہ‎ moll “s, TN 


Linear programming optimization models have been used extensively as tools 
for energy systems analysis. In these models, changes in technology and price 
structure affect patterns of resource utilization through interfuel substitution. 
However, total energy dewand and supply constraints are generally specified 
deterministically. Since price changes affect the quantity of energy used and 
the supplies that may be forthcoming, in addition to causing interfuel substitu- 
tíon, some mechanism for price responsiveness of supply and demand constraint 
values is required. in this paper, several methods of representing supply and 
demand elasticities, which were used in an energy system optimization model, 
are discussed. 

The model into which the supply: and demand elasticities were introduced is 
a linear programming model of the United States energy system (5) The problem 
is to minimize the annualized cost of satisfying the energy requirements of the 
system (which arise from demands for space heat, air conditioning, electricity, 
automobile transportation, etc.) for a particular planning year subject to con- 
straints on environmental impacts, the availability of supplies (resources and 
electric generating capacity), and the feasible range of intefuel substituta- 
bility. Time dependent features of substitutions are handled by additional con- 
straints. In some respects the problem resembles the classical transportation 
problem since what must be determined is the optimal allocation of supplies to 
demands. The representation is modified by the inclusion of efficiency coeffi- 
cients in the supply and demand constraints, environmental impact restrictions, 
and equations reflecting technical features of the energy system. Figure 1 shows 
a graphical representation of the model, and defines terms. The s supply nodes 
and d demand nodes cen include all supply and demand categories. 

To provide a feasible path between a supply and demand category, both a 


supply and utilizing technology must be identified. For a given path, 1, a 


resource S, is converted to intermediate energy form, X at an efficiency, 55 
In turn, the intermediate energy form is used to satisfy demand D, at an 
efficiency es ል cost P; and set of environmental effects, included in f 


are also defíned per unit of intermediate energy form. 


The mathematical formulation of the model is as follows: 


Minimize 


subject to 


The summation in each constraíntequation is over only those indices j that 
represent admissible intermediate energy forms for the supply category, demand 


category, or other constraint. The full linear program array thus has + 


constraints and these may be represented as: 


n 
7 A. x. € b 


ij j 1 
T 
where the A include e ,, d ,, and f , in the previous formulation and the b 
1j uj vj wj 1 
correspond to the Sy? D and B ° The t equations that augment the supply and 
demand constraints include the environmental constraints and equations which 
reflect technical features of the energy system. L is the set of less than 
or equal to constraints, E is the set of equality constraints. For a complete 


discussion of techniques for solving this problem, see Dantzig (2) and Wagner (7) 


Several methods of representing elasticities for this model were explored. 


de 


One method used to represent supply elasticities was to allow the system to 
"buy" more of a desired supply at an incremental cost, which might, for example, 
represent increased extraction costs as a resource becume more difficult to re- 


cover. The supply curve is shown in Figure 2a. The CC represent the cost of 


the resource, and it ís assumed that the cost C, is incorporated in the cost 


coefficients Py: Then the linear programming model may be reformulated as: 


min Y px, + (C, - ۷ 


j 1^"k 


0 > 7 Ax, > SUPPLY = Q + ۷ 


oem 5244178 


k 


where Ci < ር > ۰۰۰ and Q, < Q, > ,.. (only equations relevant to the elasticity 
representation have heen included. As many supply curves as desired may be 
modeled in this way although orly one is shown here). From the dual of the 
problem, the following equations are obtained (see Danzig ^) for a discussion 
of duality theory): 

-U, +ህ 


k > ٭‎ (@, ° 7 k = 1,2,... u 


where U, is the marginal value for the SUPPLY equation, and Uk is the marginal 


value corresponding to the equation limiting Wee In an optimal soluti.n, ۷ 


implies We = 41 - ኒ› since Cr - C, < CH - የ1- Therefore, by complementary slack- 


ness, when 0 < SUPPLY <Q,» U, = 0; when Q, < SUPPLY «0,1: ولا‎ = Ck ” C4; 


SUPPLY = Q ° وا‎ > የከ - “1 The corresponding quantity-marginal value curve is 


when 


shown in Figure 2b. 


The same representation can be used for demand curves. The formulation is: 


min “ P4*4 0.77 CWE 


Te = DEMAND = Qi * F Wy 


s w s ٭ ° وی‎ 


basic 


where C, > C, <... and Q < Q, <... The corresponding dual equations are: 


- > o 
"4 + U C 


and the resuit is that when q, = DEMAND, U, < Ci; when Qk < DEMAND < Q 


k+1? 
پا‎ ፦ Cy when DEMAND = Q ° U < Cc. The corresponding (DEMAND, demand marginal 


۱ 
value) curve is shown in Figure 2c. The Ck may be negative, or may be both 


positive and negative, representing for example consumer utility functions with 


(1) 


regard to demand reduction or increase from a reference level, 


Bergendahl 


6 
(6) have used this approach to model the elasticity of peak electric 


and Manne 


demand as influenced by surcharges for using electricity during these ሮ፥ . 1068, 


The step function representation above may be reformulated more compactly 


as follows, shown for the supply curve: 


min È p X +EP y 


; 3 k k 
SUPPLY = 5 Q V, 


2 7 > 1 
Ve 2 0 


where P; = 0, ۶ = Pk + (Q1 - ዒን SE - የ1)» The proof that this problem 


is equivalent to the one posed before is presented in Appendix I. 


This formlation was used by Duloy and Norton (3) to approximate a linear 


demand function, the objective being to find the equilibrium point where price 


equals marginal cost. We decided to try to find the equilibrium point where 


the amount demanded and the marginal cost would match a piecewise linear demand 


curve such as the one shown in Figure 3a. The marginal cost is defined here 38 


the negative of the marginal value, The demand marginal value curve 18 shown 


in Figure 3b. 











dems: 
functi 


solvir 


const 


where 








he basic idea behind this was that if the marginal cost of meeting 
ain demand rises , there will be pressure on consumers through price 
isms to moderate their demands. As the level of demand drops, the 
al cost decreases and the price pressures which have acted to depress 
| are relaxed. Demand may even be stimulated, perhaps by advertising, 
:s level may rise. At the equilibrium point, the forces which decrease 
is are balanced by the forces which increase them. When the demand 
‚on is piecewise linear, this equilibrium point can be found by 

ig the following non-linear programming problem. 

Consider this problem with quadratic objective function and linear 
raints: 


በሀጅ +Y (C.-C) WI Qu - Q) 


jj ke] 1 k 
r 
- Eel “kritik 
E < b, tek 
a b, i € E-(d] (1) 
E Aa - DEMAND-Q,, y - `w 
X > 0 = و1‎ n 
۱ ) 
9 > ۷ > :سیا‎ " Y باب تنب‎ 


| C1 « C, <... and Q, < Q, €<. 
One one demand function is modeled for notational convenience; the 
sentation for additional functions is the same. If (1) is rewritten 
maximization problem, the objective function may be thought of as 
izing the difference between the utility of the demands and the cost 


eting all the demands. For any particular demand, its utility curve 


mposed of segments of parabolas, each parabo'a and its successor having 


oh. 





a unique common point. The utility of two or more demands is the sum of their 












utility functions. The Kuhn Tucker conditions (5) for this problem are: 








2 ^1 Ri + Xi = b, š € L 
W +4 
k k+r = 4 - Qk k= 1, r 
ZA, x, - b, 1 € E - ۱ 
TA + y = 
"e ሜ Q. 











(Cy “OW (Q4 Y) Zk t Uq + V, = ር 


d ሂ k+1 k = 1, r 


۲ 10 + 
ER, Um * Ziemi + 3,2, + نیت‎ * 


r+k 





























where U, is the marginal value of the demand, Y is the marginal value corresponding 
to constraint limiting Wee 

From these conditions it can be seen that if a solution to the problem 
exists, then the quantity-marginal value point obtained will be on the piecewise 
linear demand curve shown in Figure 3b. The proof is presented in Appendix II. 
Since the quadratic objective function for the minimization problem is convex 
and may be written as ply + Y cy where C is symmetric and positive semi-definite, 


there are several methods for solving the problem; for a discussion see Kúñzi (5) 


Instead, we attempted to find solutions using miner modifications to the simplex 


method. 


















The basic idea was to start with a given demand vecter q^ find the 

marginal values y° from the linear programming probiem (where ba - 94 ) 
5: ኣ 

min Py j 


` - E 
5 < 0, 


and if the points (ag? ሀ4) did not lie on the demand curves, to change 4 and 


try again. If for each d, qa could be changed to the quantity corresponding 


to Ua on the demand curve without changing the basis, the solution would have 


been obtained and iteration would stop. (It is easy to see that the solutior. 
to the quadratic programming problem can then be reconstructed, since an 

optimal solution to (2) satisfies most of the Kuhn Tucker conditions, and one 
can solve for the W,, V, 889 Z, to satísfy the rest.) A natural approach is 


to change a^, the initial demend vector, to a, the vector of demand levels 
/ 


corresponding to U? on the demand curves. In problems where 94 > 94 for all d 
implies -U, >-U, የ6ኛ all 4, that is, whenever one or more demand levels increase, 
none of the demand marginal costs decreases, this method yields a sequence which 
either converges or oscillates between two limit points. Let q^ be the demand 
vector where demands are at the lowest permissible levels on the demand curve, 
Let " atis be the demand vector 


U? the associated demand marginal value vector. 


corresponding to ۳ Then a < e. 80 .v! > J^, and since the demand curves 


2 O 
are decreasing in demand, q^ < a; q^ <q by choice of q. q >q implies 
1 0 2 ها‎ . 1 ጋ 


> 


2 4 ۹ 5 ° 
-U > -Uf hence q <q. In general, q <q <q >... andq 24 q >... 


The condition that marginal costs of any demand is non-decreasing as 


demand levels rise may not hold in practical problems. For example, we found 











th 
at the marginal cost of meeting air-conditioning demand may decline as space 


heating demand increases. 


air 


generating capacity will be needed to satisfy it 


heat demand results in a winter peak, the 


When the summer-winter electric peak demand is balanced 


-conditioning demand has a relatively high marginal cost since additional 


However, if increased space 


availability of seasonal off-peak 


capacity reduces this marginal cost. 


We also found in several test problems which met the above condition 


that oscillation vetween two demand vectors occurred 
a nodifica-ion of the algorithm above, 


algorithm may not yield an exact answer to (1), however 


in the data on demand curves, 


1. 








Therefore, we adopted 


with success in all cases so far. The 


given the uncertainties 


we did not feel that grea er accuracy was necessary. 


Solve the linear programming problem 


min ۶ 7 x 





with an initial vector q^ 


At iteration k, having found U 


the demand curve to which it corresponds 









for the demands (a, = ba 4 ረ D). 


d? d € D, determine the quantíty Pi on 


k ድ 
, and set 94 ° 9, if this can 


be done without changing the basis. 


k*, k 
If 94 'ዓ4 - l| > 0.02 for de D, stop. Otherwise find the first d 


such that this does not hold, set DU 


k+1 k 
94 = 94 for the rest. 


go back to 2, 


ሂት k 
- (q; + q I) /2 for this d and 
d 


Solve the new linear programming problem and 














So far each problem has converged quickly, within 
testing purposes, the elasticities assumed were ~0.3 ! 
reference point (Q,; Po» other points on the curve we 


10/3 
P = P. (q /Q) 


for Q = 0.8 Q ° 0.9 .ی۹‎ Q ° 1.1 ዓባ 
the supply curves ranged from 0 to 1.5 9, where q. wa: 
the supply; the Cr were chosen arbitrarily. 

In the preceding, we have shown how elasticities 
be incorporated into energy system linear programming 
through more realistic modeling of the supply and dem 


yield better insights into the effects of technologic 


decisions on the energy system. 





a few iterations. For 
or each demand. Given a 
re calculated by 

The step functions for 


; a reference level for 


of supply and demand may 
models. We believe that, 
and levels, the results wiil 


al changes and policy 
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Figure 1 Graphical Representation of Linear Programming Modei 






Figure 2 Step Function Supply and Demand Curves 











Figure 3 Piecewise Linear Demand Curves 
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Figure 2c 
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Uy = 71 ۳ Q, 72 ኞ CCP say ” P )/ (01 1 - E QP, ^ 0, which in turn implies The Kuhn Tucker conditions imply that 1፻ a solution to the problem (1) 





> Vk = 1. By Lemma 1, only v, and Viel are basic, so q, v, + Ar Viel 


Q, > Q1 and the constraint SUPPLY = . Q. V, is violated. 





Iheorem: If q, < Q = SUPPLY < q, then V = (Q - Q)/( exists, then the demand-marginal value point obtained lies on the piecewise 
1 P p+1 Q 







linear demand curve in Figure 3B. 




















and p41 079)/Q, - ዒን, p> 1. 


Lemma 1: No two non-adjacent variables V, and Mey? j - 1, are positive. 


Lemma: If We > 0 then Wiel - Qe? - 41° 


Proof: Suppose not. The Kuhn Tucker conditions require that 





Lemma 3: If q, > Q > +1? two of the ve ate positive, 





Proof: Lemma 1 shows that no more than two are positive. Suppose only one 
Proof: If they were, the dual equations corresponding to y 




















and V V, is positive, Then V, = Q/Q, < 1 in order to satisfy Q V, = Qand V, > ۰ 
r+ , k y n , - - -. = 
" ! : k k k (C1 e) W Q1 ባኔ) “x d "4 7 Ye Ct] 
since < QS inequality holds. - zm. 
q, Q دوب‎ neq y holds Again, ዒ U, + 7 Pe Since 






Vk < 1 and noother v, are posítive, U, = 0 by complementary slackness. Hence / (Q ) 5 +U,+V = Ç 


ke2 ° Qi? ° "ኤ1. Pa * Vua ገ+2 








k+2 kt i k+1 





and complementary slackness ipis 9 * P/Q. But -Q..1 ሀ. + U, 2 “Br which implies ሀ= Pr %-1 < P/Q 









In addition since W, > 0, Z 


k 0. 


= 0; since W < Q - 04. , ۷ ። 
7 ኗ+1 k+1 k k+l 
ولا‎ * Pij - ESO و‎ Ys U, = =P, + çU ۰ 


Ar U, + > -P 


It follows ኩሃ Lemmas 1, 2 and 3 and the requirements that ፲: % Vk = Q, 







Bu But then U, = C - (ር 
ኒ 





i Vy <1, that if ዓ > > Sr then ^ and V 1 are basic. This means that 


- - QU) ` ۷ > 6 and 


kel Ep? v Q1 k+1 





ke]? from the dual equation corresponding to 41° 











U. > 0, as was seen in the proof of Lemma 2, hence V + Vel = ]. The unique 
therefore , l ሀ ጻር - (ር -C,.,)W /(Q - 0 ( + 7 - ር ° 
solution to the equations V + V = 1 and q, V + V *QisV = ë k+2 ws .. w 0 k ۹ -— 
-Q u, Y, + (=P, + Q U_) 2 -P 0:۴ p * 1ب‎ ۶ P 
k+l s k k 8 k+1 


(Q1 ከ 9ን/ (ዒ.+1 ° q), 1ہم۷‎ = )Q - qq ° q). 












The possibilities may now be analyzed case by case. If W = 0, then 


C ^ P! Q4 - Q.)‏ * ولا 


But ولا‎ = (P. F/G, - Q) > ወ 







DEMAND = Q+] and -Z + U = 1" oF U - ር. ° 
k+l ` P/Q, ` Qg)» due to the convexity 





of the function £(Q), where ... Wi ۳ eel ዘ Q. — ° +1 aye پت‎ ٤ I (CH ~ C) Ve Q1 " q.) ۱ 












£(Q) = 0 
E(Q) = Py + (Pq ` PD (Q 9,)/(9 






952950. ር + (ር - C_) (DEMAND - 
1+1 ~ 9%) Q 5Q > ۹ 


- Q ( which is the eruation of the line 
r+1 r+1 r 


051 / Qu 





1*1 







segment lying between (Q,> Cc) and (Qari Cup: If 7 = Q1 - Q U, = 


° ۹ € ۲۰ In general, if “> 0 then 


f i ° 
(Q) ís convex because the slopes Pr P/Q, - 91) - CH - C, are (strictly) 






increasing. 















Lemma 2: If q, < Q< Q. no two variables Vi» 1 > p*l, are positive. 


U 2 C M - - . e 
4 keg tf Wk t Qu - Qo Ug < Ck; 1۶ 0 < پل‎ < Q. - Q, then by the 
Proof: By Lemma 1, only adjacent such variables need be looked at. Again the 


dual equations and complementary slack i = ዓ- - 
mp y ckness imply ሀ (P41 P/Q) Q ) and 




















lemma, DEMAND = Q 


(C 


which is the equation of the line segment lying between (Qs C! and 


(Q1? "ቬ+1 


Thus the demand-marginal value point lies on the piecewise linear demand 





curve. 


elg "k Cea 
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r 6,,”9.) # =Q, e وا‎ $^ 
A ™ 7 5 


- Q.) = የኤ1 t (የኤጊ1 ` Ch) (DEMAND - ባ,,1)/(0ኬኤ1 - ۰ 
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